AD-hlO  457  WISCONSIN  UNIV-MADISON  MATHEMATICS  RESEARCH  CENTER  F/6  20/4 

HOMOGENEOUS  AND  NON-HOMOGENEOUS  BOUNDARY  VALUE  PROBLEMS  FOR  FIR— ETC  (U) 
NOV  81  H  B  DA  VEIGA  DAAG29-80-C-0041 

UNCLASSIFIED  MRC-TSR-2300  NL 


Cl®  FILE  COPY 


Mathematics  Research  Center 
University  of  Wisconsin-Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 

November  1981 


(Received  September  28,  1981) 


/ 

(  ; 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  O.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


Approval!  for  public  rolooso 
Distribution  unlimited 


-  1 

-  2  02  03  089 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

HOMOGENEOUS  AND  NON-HOMOGENEOUS  BOUNDARY  VALUE  PROBLEMS  FOR  FIRST  ORDER 
LINEAR  HYPERBOLIC  SYSTEMS  ARISING  IN  FLUID-MECHANICS  (Part  I) 

H.  Beirao  da  Veiga* 

Technical  Summary  Report  #2300 
November  1981 

/Xv*.  ,  i  -r  ABSTRACT 

-  .1-  i  i* 

b 

^J4e"prove  the  existence  and  the  uniqueness  of  classical  and  strong 
solutions  for  a  class  of  non-homogeneous  boundary  value  problems  for  first 
order  linear  hyperbolic  systems  arising  from  the  dynamics  of  compressible  non- 
viscous  fluids.  The  method  provides  the  existence  of  classical  solutions 
without  resorting  to  strong  or  weak  solutions.  A  necessary  and  sufficient 
condition  for  the  existence  of  solutions  for  the  non-homogeneous  problem  is 
proved.  It  consists  of  an  explicit  relationship  between  the  boundary  values 
of  u  and  those  of  the  data  f.  Strong  solutions  are  obtained  without  this 
supplementary  assumption.  ^jSee  Theorems  3.1,  4.1,  4.2,  4.3  and  Corollary  4.4? 
see  also  Remarks  2.1  and  2.4. 

In  this  paper  we  consider  equation  (3.1)  below.  In  the  forthcoming  part 
II  we  prove  similar  results  for  the  corresponding  evolution  problem. 
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SIGNIFICANCE  AND  EXPLANATION 


We  undertake  a  mathematical  study  of  the  equations  (1.1)  governing  the 
motion  of  a  compressible  non-viscous  fluid  in  some  region  of  space.  The 
equations  are  presented  in  J.  Serrin’s  article  in  the  Handbuch  der  Physik,  or 
in  the  classical  treatises  of  L.  Landau-E.  Lifchitz  or  L.  Sedov. 

A  rigorous  mathematical  study  of  these  equations  is  quite  difficult  but 
of  physical  interest  since  existence  and  uniqueness  properties  are  a  test  of 
the  validity  of  the  model  under  discussion.  Beyond  that  qualitative 
properties  of  solutions  often  underlie  the  estimates  derived  in  studying 
existence  problems.  Finally,  rigorous  numerical  approximations  can  frequently 
be  derived  from  the  mathematical  apparatus. 

A  general  approach  to  studying  non-linear  equations  is  to  use  approxima¬ 
tion  by  a  family  of  linear  ones.  In  general  different  linearizations  can  be 
found.  It  is  important  to  find  simple  ones  not  only  to  simplify  the 
mathematical  study  but  also  for  possible  use  in  applications,  particularly 
linear  numerical  approximations. 

A  study  of  the  non-linear  equations  (1.1)  was  done  in  D.  G.  Hoin's  paper 
[5]  and  in  the  author's  papers  (2),  [4].  However,  the  linearizations  utilized 
in  these  papers  are  not  the  simple  canonical  approximation  (1.2)  used  here. 

In  fact,  the  known  results  in  the  literature  are  not  applicable  to  this  last 
problem.  The  aim  of  our  paper  is  the  study  of  the  linearized  equations 
(1.2).  Actually,  we  study  a  class  of  boundary  value  problems  for  linear 
hyperbolic  symmetric  systems  which  contains,  as  a  particular  case,  the  equa¬ 
tions  under  consideration.  To  propound  the  central  rule  of  linear  hyperbolic 
systems  in  pure  and  applied  mathematic*  s«»ms  superfluous.  We  only  point  out 
that  in  the  course  of  our  study  many  i>  >ing  problems  arose  -  for 

instance,  the  compatibility  conditions  conk.  , ting  the  boundary  values  of  the 
unknown  vector  field  u  to  those  of  the  external  forces  f  (see  Section  4). 

Finally,  we  note  that  the  proofs  given  here  are  quite  simple,  and  use 
only  basic  results  in  functional  analysis  and  partial  differential  equations. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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t.  Introduction.  The  motivation  for  the  present  work  was  the  equations  of  motion  of 


.3  = 


(x2,Xj).  The  governing  non-linear  equations  are  then 


+  (v*V)v  +  p_1Vp  »  g  , 


(1.1) 


-1 


+  v*Vp  +  p(-j£)  div  v 
3S 

It 


+  v*VS  *  0  , 


( t . 0 i x ' )  -  0  , 

v( 0, x )  «  v0(x)  , 
p(0,x)  -  P0(x)  , 

S(0,x)  »  SQ(x)  , 

where  the  velocity  v  «■  *vi'v2'v3^'  the  pressure  p,  and  the  entropy  S  are  unknow 
in  ]-T,T[  x  R^.  The  functions  g(t,x),  vQ(x),  pQ(x)  and  S0(x)  are  given.  Furthermore 
the  equation  of  state  of  the  medium  p  =  p(p,S)  is  a  known  function  of  p,S  verifying 
p  >  0  and  3p/3p  >  0. 

The  motion  of  compressible  non-viscous  fluids  was  studied  by  Shin  for  small  initial 


data  (5J  and  by  us  tor  arbitrary  initial  data  [21,  13],  14] 


(1) 


However,  the 


*Oepartment  of  Mathematics,  University  of  Trento  (Italy). 

^'*Later  in  a  paper  independent  from  ours  Agemi  extended  the  approach  of  Qjin 
to  arbitrary  initial  data  (see  (1)). 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


linearizations  used  in  these  papers  are  not  the  simple  ones  which  consists  in  studying  the 
first  order  hyperbolic  system  (the  linearization  procedure  decouples  the  variable  S) 

+  (wV)v  +  aVp  «  g  , 

+  w*Vp  +  b  div  v  =  t  , 

(1.2)  v,  ( t,0,x 1 )  >  0  , 

V  (  0,  X  )  =  vQ(x)  , 

P(0,x)  «  p0(x)  . 

Here  a(t,x)  and  h(t,x)  are  given  positive  functions  and  the  given  vector  field 
w(t,x)  verifias  Wj(t,0,x')  »  0. 

As  pointed  out  in  Bain's  paper,  the  known  results  for  linear  hyperbolic  systems  do  not 
seem  suitable  for  these  problems.  For  this  reason  we  will  study  the  problem  of  the 
existence  of  classical  (i.e.  differentiable)  solutions  for  a  class  of  boundary  value 
problems  which  contains  as  a  particular  case  the  system  (1.2)  and  other  systems  arising  in 
fluid-mechanics.  This  will  be  done  in  this  paper  (stationary  case)  and  in  a  following  one 
(part  II.  evolution  case).  The  method  used  here  also  provides  a  simpler  proof  when  the 
problem  is  posed  in  the  whole  space.  In  this  last  case  the  statements  and  proofs  are 
obtained  by  dropping  all  the  assumptions  concerning  the  boundary. 

In  this  paper  we  avoid  the  use  of  weak  solutions,  mollifiers,  and  negative  norms  by 
proving  directly  the  existence  of  classical  solutions  which  are  our  main  concern,  in  view 
of  problem  (1.1).  The  existence  of  strong  solutions  then  follows  as  a  consequence.  Recall 
that  in  the  fundamental  works  of  K.  O.  Friedrichs  (7)  and  P.  r).  Lax  and  R.  S.  Phillips  [8] 
one  starts  by  proving  the  existence  of  weak  solutions. 

we  prove  the  existence  of  classical  solutions  without  assuming  the  boundary  space  N 
to  be  maximal  non-positive  and  the  boundary  matrix  to  oe  of  constant  rank  on  the 
boundary.  See  Theorems  3.1,  4.1,  4.2,  4.3  and  Corollary  4.4 »  see  also  Remarks  2.1  and  2.4. 

An  essential  tool  in  our  method  will  be  the  introduction  of  a  space  Z  of  regular 
functions  verifying  not  only  the  assigned  boundary  conditions  but  also  some  “complementary 
boundary  conditions"  such  that:  (l)  the  boundary  integrals  in  (2.10)  vanish  for  every 
u,v  «  Z;  (ii)  Z  is  dense  in  Y  (roughly  speaking:  the  complementary  conditions  loose 
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sense  in  Y)>  (iii)  there  exists  an  homeomorphism  D  from  Z  into  x  for  which  (2.28) 
holds;  (iv)  Lu  t  y,  Vu  S  Z.  These  conditions  could  be  weakened,  but  the  above  form  is 
sufficient  for  our  purposes. 

In  order  to  simplify  the  exposition  we  treat  the  problem  in  the  half  space  »?. 
However,  by  standard  methods  one  can  adapt  the  results  to  open  regular  subsets  tt  of  R™. 
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Remark  2.1.  The  results  ana  proofs  stated  in  this  paper  can  be  easily  adapted  it  the 
assumption  "for  every  x  S  rP,  H(x)  is  diagonal"  is  replaced  by  the  more  general 
assumption  "for  every  x  S  rP,  H(x)  is  symmetric,  moreover  for  every  x  e  R P_1  it  has  the 
form 


(2.4) 


H(x) 


Hp(x) 


0  “n-p'^J 

with  Hp(x)  and  Hn_p(x)  matrices  of  type  p  *  p  and  (n  -  p)  x  (n  -  p)  respectively" 
(definition  of  p  in  (2.13)).  The  assumption  (2.1)  is  then  replaced  by  "H  is  uniformly 
positive  definite".  In  this  more  general  case  the  scalar  products  (u,v)h  and  ((u,v))h 

B  #v  n 

(definition  below)  are  replaced  by  (Hu,v)  and  (Hu,v)  l  (h  4^-,  -y—)  respectively. 

J-1  (3|XJ  XJ 

Moreover  the  operator  D  in  <2. 27)  becomes  Dv  =  Hv  -  div(HVv)  '  where 
v  «  (v1,...,vn).  Note  that  now  equation  Dv  «  f  is  an  elliptic  system  of  n  equations 
instead  of  n  single  (decoupled)  elliptic  equations.  Thanks  to  the  boundary  assumption 
(2.4)  the  operator  D  is  again  an  homeomorphism  of  Z  onto  x,  moreover  (2.28) 
holds"'.  The  reader  easily  verifies  that  our  proofs  hold  again. 


Now  let  l  be  the  partial  differential  operator 


(2.5) 


Lu  =  £ 


J-1 


H-V  * 


where  u  »  (u1,...,uj))  is  a  vector  function  defined  in  rP.  In  equations  like  (2.5)  in 
which  matrices  act  on  vectors  these  last  are  always  to  be  considered  as  column  vectors. 
Note  that  equations  (1.2)  can  be  written  in  the  form 


(2.6) 


3u  v  „-1  J  3u 

li+  i  H  *  IT 

j-i  j 


where  u  »  (v, ,  v2,  Vj.p) ,  f  =  (<}  2,<3y  *■)  and 


(3) 


This  means  that  (Dv), 


(4) 


7~r hk,v. .  f 

i-1  kl  1  ,7= 


f  J 


p  3vi 

3X  (  *  hki  ajT- ' '  for  k  =  1 - - - 

T=1  J  i=l  K1  dxj 


Clearly  with  ( (u , v) ) h  replaced  by  the  scalar  product  indicated  above. 
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0  w,  0  0 


0  0  w,  0 


lb  0  0  w, 


h-’a2  - 


Let  us  return  to  the  general  case  (2.5).  Let  us  define  AJ  =  H_,AJ  hence 

a^  ■  h”’a^k.  For  convenience  we  will  use  the  notations  lAJl  ,  =  maxlaJ.l 
,  C  i.k  1  C 

>Al  ,  s  max  I A  I  ,  and  similarly  for  the  other  matrices  used  in  this  paper.  By  definition 
C  J  C 

•  hi  .  =  maxlh.  I  ,  »  IH)  Let  be 
c‘  k  k  c*  c1 

X  5  (L2(j»"))n,  V  =  {H1(R®)ln.  Z  s  (H2(R®)Jn 
and  define  the  scalar  products  in  X 


(u.v)  =  l  j  u^dx,  (\J,  v)h  =  i  /  nkvkhk( 


k'1  R- 


k-1  Rm 


and  also  the  corresponding  norms  |*|  and  |*|h  (which  are  equivalent)  and  the  scalar 
products  in  V 

((u.v))  3  (u.v)  +  (Vu.Vv),  ((u,v))h  3  (u,v)h  +  (Vu,Vv)h 
and  corresponding  (equivalent)  norms  •*•  and  !•!  .  By  definition 

IB  3y 

(Vu.Vv)  =  £  (-y— — ,  v — )  .  To  point  out  that  X  is  endowed  with  the  norm  |»|  we 

J-1  XJ  XJ  h  h 

sometimes  write  Xh  instead  of  X.  A  similar  remark  holds  for  V. 

For  the  reader's  convenience  we  state  the  following  two  lemmas  (which  proofs  are 
classical): 


Now  let  p  be  an  integer,  0  <  p  <  n  and  let  us  define 
N  »  (u  6  R°  !  -  •••  -  Up  -  o),  N'1  -  lu  ff  bF*  :  up+1  »  •••  »  U„  »  O), 

PNu  ■  (u  t,...,un),  P^|u  *  (u^,...,up).  We  assume  that  the  boundary  conditions  are  given 

by  P^u  «  0  on  the  boundary  or  more  explicitly 
N 

(2.13)  u K(0,x')  -  0,  k  «  1,...,p,  for  x1  6  R®'1  . 


we  also  assume  that  the  boundary  matrix  A1  ■  njAJ  verifies  for  each  x'  6  R®-1  the 

J-1 

following  assumptions  (the  reader  is  also  referred  to  the  papers  [7]  and  (S);  see  in 
particular  the  sections  5  and  8  of  [7])t 


(2.14) 


a’(n)  c  n1 

A1  (N1)  C  N  . 


Moreover  we  assume  that  for  each  x  8  Rn~1 

(2.15)  AJ(N)  C  N,  J  m  2 . a  . 

With  the  boundary  conditions  in  canonical  form  (2.13)  assumptions  (2.14)  and  (2.15) 
are  respectively  equivalent  to  the  following  ones:  The  matrices  AJ  take  on  the  boundary 
the  form 


(2.16) 


.1 


rRJ 


J  «  2, 


where  M(x*).  RJ(x‘)  and  SJ(x')  are  p  x  (n  -  p),  p  *  p  and  (n  -  p)  *  (n  -  p) 
matrices  respectively.  MT  denotes  the  transpose  of  M.  By  preceding  assumptions  RJ 
and  SJ  are  symmetric.  The  ranks  of  these  matrices  are  free. 

Note  that  if  p  »  0  or  if  p  ■  n  the  conditions  (2,15)  disappear  and  condition 
(2.14)  becomes  A1  a  0  on  the  boundary:  in  case  p  *  0  it  suffices  that  A1  was  negative 
semi -definite;  see  Remark  2.4.  Note  also  that  in  the  particular  case  (1.2)  one  has 
p  -  1  and  the  matrices  (2.7)  verify  the  assumptions  (2.16). 

Finally,  we  describe  the  assumptions  on  the  lower  order  term  Bu  in  equation  (3.1). 

We  assume  that  B  is  an  n  *  n  matrix  valued  function  defined  in  R®  and  verifying 

(2.17)  bik  S  c’(R“),  i,k  -  1 , . . .  ,n  . 
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Moreover  we  assume  that  for  each  x  on  the  boundary  one  has 
(2. IB)  B(N)  C  N  . 


This  last  condition  means  that  on  the  boundary 


<2.  19) 


uB2 


where  Box'll  B2(x')  and  B^tx1)  are  matrices  of  types  p  x  p,  (n  -  p)  *  p  and 
(n  -  p)  x  (n  -  p)  respectively.  Clearly  b  6  L[Y;Y ];  see  definitions  (2.22). 
we  now  define 


(2.20) 


max (sup 
uCX 


|  (Bu,u)h  +  -  a(u,u) | 


l» 


| ( ( Bu , u ) ) h  +  j  B(u,u ) | 

sup  - - - }  . 

u6¥  luC 

h 


From  our  assumptions  on  the  coefficients  it  follows  the  existence  of  a  constant  c 
such  that 


(2.21  ) 


I A I  . Ihl 


Moreover  we  define  the  following  spaces: 

Y  5  {u  G  y  :  u  e  N  on  the  boundary}  -  IhJ(R?|)]p  x  [H'(K™))n“P  , 

and 

Z  z  {u  G  Z  :  u  G  N  and  6  on  the  boundary}  »  x  (H^(H®))n-P  . 


Note  that  Z  is  dense  in  Y  and  Y  is  dense  in  X. 

Before  going  on  we  discuss  in  the  following  remark  some  features  of  the  boundary 
assumptions. 

Remark  2.4.  Assume  that,  instead  of  (2.1b) )f  the  boundary  matrix  has  the  more  general 

form 

0  m" 

T 

M1  Q 

where  y  is  a  symmetric  and  negative  semi-definite  matrix;  this  is  equivalent  to  replace 
condition  (2.14)1  oy  condition 
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(2.22) 


Vu  6  N 


vx  •  e  r®“' 


u'a'u  <  0, 

Linda r  this  weaker  assumption  the  existence  results  proved  in  our  paper  (including  those  of 
Section  4)  holds  again  provided  that  one  assume  that  A  >  A^  instead  of  |A|  >  Aq.  The 
proofs  remain  unchanged  in  the  essential  features. 

(ii)  Energy  estimates  are  obtained  under  weaker  assumptions  than  those  needed  to  get 
existence.  The  a  priori  bound  |u|),  *  (|Aj  -  A^)  1  |  Au  +  Lu  +  Bu|,,,  Vu  a  Y,  holds  without 
assumptions  (2.14>2,  (2.15)  and  (2.19);  moreover  assumption  ( 2 . 1 4 ) 1  can  be  replaced  by 

(2.22)  if  we  take  in  account  only  the  values  A  for  which  A  >  A^.  Analogously  the  a 
priori  bound  *u*h  *  (|A|  "  A0)_1IAu  *■  Lu  +  Bul^,  Vu  6  Z,  holds  without  the  assumptions 
(2.15)  and  (2.19);  moreover  if  we  take  in  account  only  values  A  >  AQ,  the  assumptions 
(2.14),  and  (2. 1 4 )  2  can  be  replaced  respectively  by  (2.22)  and  by  u*a’u  <  0,  Vu  6  N*, 

Vx '  S  R®"1.  The  remaining  assumptions  (2.14)2,  (2.15)  and  (2.19)  are  utilized  only  to 
get  L(Z)  C  Y  and  B(Y)  C  Y;  it  seems  clear  that  these  assumptions  could  be  weakened  if 
one  only  wants  to  prove  existence  for  strong  solutions. 

Lemma  2.5.  Assume  that  (2.14),  holds .  Then 

(2.23)  (Lu,v)h  -  -(u,Lv)h  +  a(u, v),  Vu,v  6  Y  . 

In  particular 

(2.24)  |  (Lu  +  Bu,u)h|  4  \>lulh'  Vu  C  Y  . 

In  fact,  under  the  hypothesis  of  the  lemma  the  last  term  in  equation  (2.8)  vanishes. 

Assume  now  both  hypotheses  (2.14)  and  let  u,v  C  Z.  If  t  *  1  the  tangential 

derivatives  4~  and  x-~—  belong  to  N  and  from  (2.14),  it  follows  that  the 
3xt  3xt 

corresponding  integrals  in  equation  (2.10)  vanish.  This  also  holds  for  t  »  1  as  one 
shows  by  using  (2. 14)2  and  by  recalling  that  the  normal  derivatives  and  belong 

to  N*.  Hence  one  gets  the  following  result; 

Lemma  2.6.  Assume  that  (2.4)  holds.  Then 

(2.25)  ( (Lu, v) )h  «  - ( ( U , Lv ) )  +  B(u, v),  Vu,V  6  Z  . 

In  particular 

(2.26)  | ( (Lu  +  Bu,u))h|  <  AQlul3,  VU  B  Z  . 
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3.  Classical  Solutions  for  f  6  y.  Strong  Solutions.  The  assumptions  in  this 
section  are  (2.1),  (2.2),  (2.3),  (2.14),  (2. IS),  (2.17)  and  (2.18).  Recall  also  the 
definitions  (2.S)  and  (2.20).  The  boundary  conditions  are  given  by  (2.13).  Under  these 
assumptions  we  prove  the  following  theorem  of  existence  and  unicity  for  classical  and 
strong  solutions  of  equation  (3.1)^: 

Thereom  3.1.  Let  the  above  conditions  hold  and  let  X  6  R  be  such  that  I X I  >  X  . 

- - -  I  I  Q 

Then  (i)  for  each  f  C  Y  the  equation 

(3.1)  Xu  +  Lu  +  Bu  »  f 

has  a  unique  classical  solution  u  8  Y.  Moreover 

(3*2)  ,u,h  4  w±-r0  ,f,h  • 


(ii)  For  each  f  6  X  the  equation 
(3.3)  Xu  +  Lu  +  Bu  =  f 

has  a  unique  solution  u  6  D(L).  Moreover 

<3,4)  iuih 4  ixr^-x-  lflh  * 

Proof.  We  give  two  different  approximations.  By  Galerkin's  method  and  by  elliptic 
regularization. 

1st  method.  I^t  l*a)«  *  8  N,  be  a  base  for  Z  and  put 

(3.5)  u(t)  -  f  Vi  6  W  , 

8*1 

Select  the  real  numbers  aa  the  solutions  of  the  linear  non-homogeneous  system 

(3.6)  ((Xu(t),a  ))  +  <<(L  ♦  B)u(t),a  )  >.  »  (<f,a_>),  ft  <  (  . 

r  n  r  n  r 


(5) 


See  also  Remark  2.4,  (i) 


for  each  t  6 


this  problem  is  uniquely  solvable.  By  using  (2.26)  and  (2.20)  it 


2th  method.  Let  £  >  0  be  a  parameter  and  look  for  u£  tf  Z  such  that 

(3.10)  M(uC,v))h  +  (((L  +  B)uE, v) )  ±  e((uE,v))z  m  ((t.v))h,  Vv  S  Z  . 

Me  take  +e  if  X  >  XQ  and  -£  if  X  <  -XQ.  The  left  hand  side  of  the  above 

equation  is  a  bilinear  continuous  and  coercive  form  over  Z,  consequently  by  a  classical 

e 

result  of  P.  D.  lax  and  A.  N.  Milqram  there  exists  a  unique  solution  u  of  problem 

(3.10) .  By  taking  v  •  uE  it  follows  that  (  I X I  -  X  )luEl^  +■  eluElf  <  If  I  luEl..  Hence 

on  l  n  n 

(3.11)  '“e,h  <  ,f,b  ' 

G  2 

and  from  the  first  estimate  above  it  follows  that  elu  I  is  bounded  by  a  constant 

independent  of  £.  Thus  there  exist  a  subsequence  uE  and  an  element  u  6  Y  such  that 

£ 

u  — m  u  weakly  in  Y  , 

(3.12) 

elu  *z  *  0  > 
e 

when  e  ♦  0.  In  particular  (L  +  B)u  — x  (L  +  B)u  weakly  in  Y.  By  using  now  (2.28)  we 

£  G  £ 

write  (3.10)  in  the  form  (Xu  +  (L  +  B)u  -  f,l)v)  »  *  e((u  ,v))  ,  Vv  «  Z.  By  passing  to 

z 

the  limit  when  e  ♦  o  it  follows  that  Xu  ♦  (L  +  B)u  -  f  »  0.  The  remainder  of  the  proof 
is  as  in  the  1st  method  above.  Cl 
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-r~ ,  ...  m 

.V,- 


4. _ The  Non-Homogeneous  Problem.  For  convenience  define 


Hp  ^  tH'd"-1))? 


for  e  »  or  and  p  S  N.  In  this  section  we  consider  non-homogeneous  boundary 
condition  P^u  »  w  or  sore  explicitly 

(4.1)  VO,*')  ■«k(»'l1  k  -  1....,p,  for  x*  «  B**-1  . 

Let  f  6  X  and  w  6  H^/2  .  we  said  that  u  6  X  is  a  strong  solution  of  problem 

(3.1) .  (4.1)  if  there  exist  sequences  u^  6  ^ 1  *  6  X,  te  U,  such  that 


xu“>  ♦  *<*>  -f(t),  p yi) 


w  on  the  boundary  (in  the  usual  trace  sense)  and 


■  “»  ♦  u,  f***  ♦  f  strongly  in  X. 

Theores  4. 1 .  Let  f  6  X  and  w  S  H^2  be  given,  and  let  |l|  >  XQ.  Then  there 
exists  a  unique  strong  solution  u  of  problem  ( 3.1 ),  (4.1).  Moreover 


Mh  <  TX-pV  11*1  +  -  >o  *  o  *  ,b'po,,w,hi/2} 


Proof.  Consider  a  linear  continuous  operator  w  +  w  f roe  H^/2  into  y  such  that 
P^w  «  w  on  the  boundary  B*-1  (in  the  trace  sense)  and  PNw  3  o  in  t*.  By  carrying  out 
the  change  of  variables  u  -  w  ♦  v  and  by  using  Theoree  3.1  the  result  follows  easily.  □ 
For  classical  solutions  a  corresponding  result  fails.  We  start  by  giving  a  necessary 
condition,  later  on  we  will  prove  its  sufficiency.  Define  smooth  solution  as  a  classical 
solution  which  belongs  to  Z.  Let  now  toy  and  let  u  be  a  smooth  solution  of  equation 

(3.1)  with  the  non-hoaogeneous  boundary  condition  (4.1).  At  each  point  x  6  n"*1  one  has 

f  • 


piH'v  -  hpvv 

pjB  -  B,pJ  . 


where  the  operators  act  on  coluan  vectors  of  if1  and  H^1  is  the  p  *  p  diagonal  matrix 
with  diagonal  elements  hj^(x),  k  •  1,...,p.  Hence  by  restricting  equation  (3.1)  to  the 
boundary  and  by  applying  P^  to  both  sides  one  gets 


«  Ju  3u 

Mp  . tt)  "  F 
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{ 


where  by  definition 

(4.5)  P  5  P^f  -  Fx[w) 

and  is  the  operator 


(4.6) 


F^Iw]  =  Xw  + 


-1_J  3w 


J-2 


H  R 
P 


TT  * 


B.w 


which  acts  on  vector  fields  w(x')  -  (w^ (x1 ), . .. ,wp(x ' ) )  defined  on  the  boundary  k"-1.  Hence 
(4.7)  F(x)  6  range  of  Hj^(x}M(x)  a.e.  on  the  boundary  , 

is  a  necessary  condition  for  the  existence  of  a  smooth  solution.  In  order  to  reverse 
condition  (4.7)  we  want  to  state  it  in  teres  of  the  functional  spaces  used  in  this  paper. 
One  has  the  following  results 

Theorem  4.2  (necessary  condition).  Let  the  data  f  (  )  and  w  6  H 2^2  be  given  and 
define  F  by  equation  (4.5).  If  there  exists  a  smooth  solution  of  problem  (3.1).  (4.1) 
then  the  equation 


(4*8)  H;1M(gp+1...-.gft)  -  (F, . Fp) 

admits  at  least  one  solution  g  6  H^/2. 

Proof.  Note  that  F.  6  Lltd/2iHV2].  Hence  P  S  Hl/2.  Let  now  u  C  Z  be  a  solution 

3  *u  ,  3u  P 

of  (3.1),  (4.1).  Clearly  P(I('5~_)  “  •  •  • »  U~)  ®  H^/2  on  the  boundary.  Moreover 

from  equation  (4.4)  one  shows  that  pl|('5~)  i»  *  solution  of  equation  (4.8).  □ 

We  will  now  prove  in  the  next  theores  that  the  necessary  condition  stated  in  Theorem 
4.2  is  a  sufficient  condition  in  order  to  get  a  classical  solution  u  6  y. 

Theorem  4.3  (sufficient  condition).  Let  f,  w  and  F  be  as  in  Theorem  4.2.  If 
equation  (4.8)  admits  a  solution  g  6  W’/2  then  problem  (3.1),  (4.1)  is  uniquely  solvable 
in  V  (hence  in  Y  if  w  =  0).  Moreover  the  following  estimate  holds 


clhl 


v2 


(4.9) 


,u,h  *  |X|"'-  XQ  It’  +  K(2|l|  -  *0  ♦  IA'  ,  ♦  IBI  ,)]lfl  + 


+  [  2  |X  I  -  »  ♦  1*1  +  IBI  J[l  +  K(  |  X  |  +  IAI  +  IBI  ,)]lwl  ) 

u  r»  1  ^  ■  i  *3/ 2 
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where  K  is  any  real  number  verifying 


Iql  <  KlU 

9  tfV2  MV2 

n-p  p 


Proof .  For  each  i fr  denote  by  and  Y1  #  the  values  (in  the  usual  trace 

3iii 

sense)  of  i)>  and  on  the  boundary.  It  is  well  known  that  there  exist  right  inverses 

Y*’  6  L(H3/'2(Bn",);H2(*")l  and  y"’  «  L  W  1/2(Rm-1 )  ;H2<«") )  of  Y0  and  Y, 

respectively.  Define  u  ■  Y  Vi  J  •  1,...,p,  and  u  =  Y,  J  =  p  1,...,n.  Clearly 
J  U  J  J  1  J 

lul-  <  c(lwl  +  Igl  ).  Carry  out  the  change  of  variables  u  »  v  +  u.  Equation 

1  H  H  /2 

P  P 

(3.1)  becomes 


(4.11)  Xv  +  Lv  +  Bv  «=  f1  =  f  -  (Xu  +  Lu  +  Bu)  . 

Obviously  f1  Q  y.  by  using  (2.16),  (2.19)  and  (4.8)  one  easily  gets  f^  6  N  on  the 
boundary.  Hence  f^  G  Y.  By  Theorem  3.1  equation  (4.11)  has  a  unique  solution  v  G  Y, 
moreover 

(4.12,  ,v,h  <TT|_I__  ,fi,h  . 

This  means  that  equation  (3.1)  admits  a  solution  u  6  V  verifying  (4.1)  and  verifying 
the  estimate 

lulh  <  lhl’^2(»ill  ♦  [Tp~l~  +  Du  +  Bu  -  f  l)  . 

Recalling  the  above  estimates  and  recalling  also  that 

,FXlW,,H1/2<C(lXl  +  “'.1  +  ,B*  i),W,u3/2 

H  n_ 

P  P 

and  that  IP  £•  1/?  *  c*f*#  the  estimate  (4*9)  follows  with  straightforward 

V 

calculations.  The  uniqueness  is  obvious.  □ 

Consider  now  the  matrix  M(x).  We  say  that  rank  M(x)  «  p  uniformly  tor  x  6  ll*1'1 

if  the  sum  of  the  squares  of  the  determinants  of  order  p  contained  in  M(x)  is  bounded 

below  by  a  positive  constant  independent  of  x,  i.e.  if 

(4.13)  l  M2  »  d2  >  0,  Vx  6  R1""1  . 

OS  I 
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Corollary  4.4.  Assume  that^ 8  * 

(4.14)  rank  M(x)  »  p,  uniformly  on  R°-1  . 

Than  to  each  pair  (f,w)  6  y  *  H2/2  corresponds  a  unique  classical  solution 

problem  (3.1)#  (4.1).  Moreover  u  verifies  (4.9)  with  a  value  K  such  that 

f  K  <  £_  IA«  Ihl  .  if  p  -  1  , 

I  dz  c1  c1 


(4*15) 


u  e  y 


of 


1C  <  £_  lAl2j^3IAl2  Ihl  ,  if  p  >  1  , 

dz  c°  c'  c1 


where  the  norms  concern  i'  i-  boundary  spaces  C*(^_1),  t  •  1,2. 

Proof.  Let  i  ov  <»'  <H>t  of  all  p-tuples  of  integers  a  -  (<^,...,0  )  such  that 

1  *  a,  <  <  *  A'  • i  •  ,  Denote  by  A,,  J  «  1#...,n  -  p  the  J-colunn  of  the  matrix 

1  {  J 

M,  by  M  ,  a  -(»...,.>)  c  I,  the  value  of  the  determinant  whose  columns  are 
a  p 

A  ,...,A  \ by  *1'(P)  the  value  of  the  determinant  obtained  from  N  by  replacing 

<*„  a  >  o 

i  p 

the  J-column  A^  by  Hp(p, ,. . . ,Fp)  -  <h1F1 , . . . ,hpFp).  Let  finally  ga  £  (g, , . . . ,gn_p) 

^  i 

be  the  vector  column  such  that  ga  »  M^fF)  for  i  »  1,...,p  and  =  0  for 

k  0  {a  ,...,a  }.  From  the  definitions  one  has  Mg  -  f  A  g  and  arguing  as  in  the 

**  i»1  i  1 

proof  of  Crammer's  rule  one  shows  that  Mg_  »  M  H  (F  , . . . ,F  ).  The  vector 
*  a  a  p  1  p 


(4.16) 


(  l 

aCI 


aei 


M  g 
u^o 


verifies  equation  (4.8)  i.e.  the  vector  g(x)  is  a  solution  of  (4.8),  for  x  S  rP_1. 

Obviously  g  6  Finally  estimate  (4.15)  follows  from  (4.16)  with  straightforward 

calculations.  Recall  that  Knl  ,  ,  <  cKI  ,  .  In<  .  ,  •  □ 

H1/2(mB-1)  c’d""1)  H,/2(*B_1) 

Remark  4.5.  Condition  (4.7)  determines  the  linear  subspace  of  data  (f,w)  for  which 
a  classical  solution  u  of  problem  (3.1),  (4.1)  exists. 


(8) 


Note  that  condition  (4.14)  can't  be  verified  if  2p  >  n. 


.  eft.  " 
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For  simplicity  assume  the  homogeneous  boundary  condition  w  =  0.  By  neglecting  w 
the  above  linear  subspace  becomes* 

Y1  »  If  S  y :  P^f(x)  6  range  of  H"1(x)M(x)  a.e.  on  the  boundary}  . 

If  we  want  to  solve  (3.1)  for  every  y  6  V  (this  means  Y1  «  y>  condition  (4.1) 
becomes  rank  M(x)  =  p  on  the  boundary.  In  the  other  extreme  case,  namely  rank 
M(X)  =  0  on  the  boundary,  condition  (4.7)  says  that  (smooth)  solutions  can  not  exist  if 
f  (I  Y.  In  intermediate  cases  for  which  rank  K(x)  5  q,  0  <  q  <  p,  explicitly  necessary 
and  sufficient  conditions  could  be  obtained  from  equation  (4.4)  and  from  Theorem  4.3.  Tl 
was  specified  only  for  q  *  p  (Corollary  4.4)  because  this  is  the  situation  in  problem 
(1.2)  where  M  »  (0  0  1),  consequently  q  *  p  -  1. 


(9) 


Note  by  the  way  that  Y  C  Y1  C  y, 
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